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Abstract 

We performed electric and thermoelectric transport measurements of bilayer graphene in a mag- 
netic field up to 15 Tesla. The transverse thermoelectric conductivity a xy , determined from 
four transport coefficients, attains a peak value of a X y,peak whenever chemical potential lies in 
the center of a Landau level. The temperature dependence of a xyiPea k is dictated by the dis- 
order width W L . For k B T/W L <0.2, 

Q xy,peak is nominally linear in temperature, which gives 
a xy,poak/T = 0.19 ± 0.03nA/K 2 independent of the magnetic field, temperature and Landau Level 
index. At I^bT/Wl >0.5, a xyiP eak saturates to a value close to the predicted universal value of 
4 x (ln2)kee/h according to the theory of Girvin and Jonson. We remark that an anomaly is found 
in a xy near the charge neutral point, similar to that in single-layer graphene. 

PACS numbers: 65.80.Ck, 73.63.Bd, 72.15.Jf, 73.63.-b 
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Graphene, a single atomic layer of carbon sheet, has attracted considerable attention 



in the past few years Its unusual massless Dirac-fermion excitations, deriving solely 
rom its honeycomb structure, have already unveiled many interesting transport phenomena 
2|-|4j. More recently, there is a growing interest in the bilayer graphene partly due to the 
observation of tunable band gap by breaking inversion symmetry using external bias ^J. 
In addition, its chiral massive excitations are unique which give rise to the Berry phase 
2tt jsl and the two- fold orbital degeneracy in the zero-energy Landau level (LL). Several 
authors have pointed out that the coulomb exchange interaction within a LL with multi- 

n 

degeneracy can break the symmetry and induce a spin-polarized quantum Hall state [7j. 
This becomes more pronounced in a bilayer graphene in which its zero-energy LL contains 
eightfold degeneracy (two spin, two valley and two orbital degeneracies) [8|. The lifting of 



spin degeneracy in the zero-energy LL has been uncovered in single-layer 



9] and bilayer 

graphene [lOj providing the evidence for exchange enhanced Zeeman splitting. 

A fundamental issue in the quantum Hall effect (QHE) is the role of the edge current 
contribution to the observed quantized Hall conductivity in two dimensional electron system 
(2DES) with finite dimensions [n]. A current density J generated by an electric field E and 
a temperature gradient (— yT) can be expressed by J = erE + a(-yT), where a and a are 
the conductivity tensor and thermoelectric conductivity tensor, respectively. By assuming 
no transverse temperature gradient in the steady state, a can be expressed in terms of 
conductivity cr xx , Hall conductivity <r xy , thermopower S xx and Nernst signal S yx as 

c^xx = CxxSxx cr xy Sy X andtt X y = (XxySxx + cr xx S yx . (1) 

We have used the definition of S xx = — E x /| y T| (positive for hole) and S yx = E y /| y T|. 
Under an intense magnetic field, there exists current-carrying edge states within a magnetic 



length £b = yh/(eB) to the boundary. In the presence of E, the induced difference in the 
edge currents at opposite sides due to the shift in chemical potential /x can be calculated to 
give <5I e dge = z/ e 2 Vn/h, where v is an integer and Vh is the Hall voltage. This agrees well 
with the observed quantized Hall conductivity o"h = ^e 2 /h with v representing the filling 
factor. An alternative way to generate a non-zero 5I e dge is to apply a (— yT). The resulting 
Pledge is also quantized, giving a universal value in the transverse thermoelectric tensor 
a xy = (ln2)kRe/h whenever fi sits at the center of a LL according to Girvin and Jonson's 



(GJ) theory 12j,|l3j. Recently, measurements on single- layer graphene 14| found that a 



•xy 



is consistent with the predictions of GJ. In this work, we report electric and thermoelectric 
transport measurements in bilayer graphene. a xy is determined from cr^, <r xy , S xx and S yx . 
It reaches a peak value of a xyiP eak whenever /i crosses the center of a LL. The temperature 
and ma gnet ic field dependences of a xy pea k are then discussed and compared to the Mott 
relation [15j . 

Bilayer graphene samples were prepared by mechanical exfoliation from high quality 
bulk graphite on a doped silicon substrate, which serves as a bottom gate, with 300 nm 
of silicon dioxide. The number of layer was first identified from the optical contrast and 
then further confirmed by the pattern of a xy plateau in quantum Hall regime. The device 
shown in the inset of Fig. [2] was fabricated using electron beam lithography 
comprising four electrodes, a local heater and two local thermometers, followed by Cr(l 
nm)/Au(40 nm) film deposition and lift-off process. A (— yT) across the bilayer graphene 
was modulated by feeding an alternating current to the local heater at a frequency uj ~ 3 
Hz. The thermoelectric signals S xx and S yx are then detected at second harmonic 2u using 
a lock-in amplifier. Similarly, the temperature gradient is determined by feeding a direct 
current to the thermometer while detecting its resistance oscillation at 2u>. The (— yT) 
we used is typically in a range from 10 mK//im to 50 mK/yum. Both and S yx signals 
are found to scale linearly with (— yT) as expected. The results from two devices, labeled 
as S37 and S45 with different sample geometry, show consistent behavior which would be 
described in detail as follows. 

Fig. [1] shows Sxx as a function of applied gate voltage V g at different temperatures ranging 
from 15 K to 300 K in S37. S xx equals zero at charge neutral point (CNP) (Vg = V g — V g o = 
0, V g o=+5 V) while its magnitude increases rapidly away from CNP and reaches a maximum 
value of S m at \V' g \ ~ +15 V corresponding to a carrier density n c « 3 x 10 12 cm -2 determined 
from its Hall coefficient. S xx monotonically increases with increasing temperature giving 
|S m | 95;uV/K at 300 K. The upper-right inset of Fig. [1] plots the temperature dependence 
of |S m | and |Sxx(V g = 80V) |, in which T 1 / 2 and T, respectively, power law fittings shown 
as red lines agree well with the data points. On the other hand, 1/Sxx is practically linear 
with V g away from CNP (|V g | >20 V), which is shown at selected temperature of T=300 
K, 120 K and 50 K in the lower-left inset of Fig. [TJ From semi-classical theory for a 
free electron system, the thermopower can be described by S xx = — 7r 2 /2(kB/e)(kBT/ep), 
where e is the Fermi energy. Since ep oc k| oc V g in a bilayer graphene, the observation of 
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S xx oc T / Vg away from CNP shown in the inset of Fig. [T] is then consistent with the free 
electron model. We also note that the T 1//2 dependence of Sm near CNP is different from 



the T-linear dependence found in single-layer graphene [14], [16|. As it turns out, the failure 
of the semi-classical theory near CNP persists to the quantum Hall regime. 

As a standard practice, <7 xx (<7 xy ) and S xx (S yx ) under a field were (anti) symmetrized with 
respect to the magnetic field to exclude the contribution from mis-alignment of electrodes. At 
15 Tesla, the observed Hall plateau of <7 xy = ue 2 /h in S45 (the upper panel in Fig. [2]), where 
v = ±4, ±8 is the filling factor, indicates the onset of quantum Hall effect even at tempera- 
ture up to 100 K. The 8e 2 /h step in er xy across the CNP (V g o ~+29.4 V) further confirms the 



bilayer graphene nature in S45, where its LL energy, given by E n = hx>cy N(N — 1) with the 
LL index N and cyclotron frequency uq = eB/m*, has two- fold orbital degeneracy (N=0 and 
1) at zero-energy LL. Using m*=0.054m e [8J], the LL energy difference AE = E 2 — E ~530 
K. From low-field Hall measurement, we determined the mobility /i c ~ 2, 600cm 2 /V-sec in 
S45 and the gate capacitance C g ~110 aF//im 2 . If naively using E=(C g /e)(7rh/2m*)V g , the 
LL spacing equals AE ~560 K, obtained from cr xx peaks in Fig. [2j which is reasonably close 
to the theoretical value (530 K). Similarly, the disorder width Wl, defined as the full-width 
at half-maximum of the peaks (lower panel in Fig. [2]), gives ~180 K. Therefore, the 
criteria of AE 3> Wl, keT for quantum Hall regime is justified. We remark the appearance 
of double-peak feature in cx xx with slanted a xy plateau near CNP, which suggests the lifting 
of spin degeneracy at zero-energy LL as pointed out recently by Zhao, et al. in a bilayer 
graphene [181]. The sheet resistance at CNP R turns out to be nearly T independent below 
30 K in zero field (the middle inset of Fig. |2]). In contrast, the Ro at 15 Tesla increases 
by 50% as T drops to 8 K, which basically agrees with previous reports of field-induced 



insulating state 



18 



19]. 



Fig. |3](a) shows V g dependence of S^ and S yx in S37 at // H=15 Tesla and T=15 K shown 
as thin black line and thick red line, respectively. The corresponding (thin black line) 
and <r xy (thick red line) vs. V g are plotted in Fig. |3] (b) with N up to 4. A double Gaussian 
fit to the (Txx near CNP shown in the thick blue line indicates an excess conductivity near 
CNP. Both Sxx and S yx exhibit repeatable oscillations with V g that effectively shifts the 
position of fi in bilayer graphene. When /i crosses the center of a LL, Sxx attains a peak 
value that progressively reduces in magnitude at higher N while S yx , on the other hand, is 
nearly zero. Around CNP (V g o — +5 V), however, a large and broad peak of S yx — +18/x 
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V/K is observed instead. Using Eq. [TJ the obtained and function of V e are 

shown in Fig. [3](c). a xy is positive and reaches a local maximum value of a xyiPea k — 3 n A/K 
for N=0(1), 2, 3 and 4 when /i crosses the center of a LL, while a xx approaches zero instead. 
We also performed a current annealing (j ~ 3 x 10 8 A/cm 2 ) on S37 , which causes the shift of 
V g o from +5 V to - 13 V as demonstrated by dashed lines in Fig. [3J Nevertheless, no other 
major variations in V g dependence are found to result from the current annealing process. 

In order to further explore the temperature dependence of a xy , we performed measure- 
ments at different temperatures up to 100 K. At low temperatures and 15 Tesla, both a xy and 
S yx are practically linear in T as demonstrated in FigJU(a), where a xy /T-V g at T=15 K, 20 K 
and 25 K all collapse onto a single curve, and so does S yx /T-V g . According to semi-classical 
theory (generalized Mott relation), a xy and S yx can be described by a xy = 2 j^-T(^|| il ) (E=M 
and S yx = ^^fT( 8t ^ H ) e=At , where tan#n = o"x y /c>xx is the Hall angle. Therefore, we expect 
that « xy /T oc (^)(^) £=M . Similarly, S yx /T oc (^fX^W The calculated ^ and 
d^y° H are shown as orange and green dashed lines, respectively, in Fig. H] with arbitrary 
units, which qualitatively agree with the measured a xy /T and S yx /T except near the CNP. 
The term (^f) e=j[1 is proportional to the density of state that gives minor influence to the 
behavior. We noticed the apparent electron-hole asymmetry in the magnitude of the calcu- 
lated and «s w hich may result from the shorting of voltage leads to the "hot-spot" 
near the current leads |20|. On the contrary, S xx and S yx , generated by a (— VT) across 
bilayer graphene, exhibit much less asymmetry due to the hot-spot shorting effect. The 
deviation from the calculated value is evident near CNP and more dramatic in S yx . The 
temperature dependence of a xy)Pea k in S37 and S45 at 15 Tesla and N=0(1), 2, 3 and 4 are 
shown in Fig. [4] (b), where closed symbols and open symbols refer to the values in S37 and 
S45, respectively. Even though there is a minor variation in a xyiPea k at different Ns and 
different samples, it basically follows T linear dependence below 30 K (k B T/WL —0.2). The 
corresponding values of a xyiPea k/T are plotted as a function of /ioH in the inset of Fig. |4] (b), 
averaged to a constant value of a xyiPea k/T = 0.19 ± 0.03nA/K 2 independent of /^oH, N, and 
T. As T increases above 80 K (IcbT/Wl — 0.5), it flattens up to a nearly constant value of 
« xy , P eak — 8.5±1 nA/K, close to the universal value of 4 x (ln2)kBe/h (red thick line in Fig. 

H. 

The disorder width Wl, estimated to be ~150 K and ~180 K in S37 and S45, respectively, 
turns out to be an important parameter for the behavior of a xy shown in Fig. H] (b). At low 
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temperature keT <C Wl, the number of disorder-induced extended states participating the 
diffusive transport increases with T, which gives rise to the T linear dependence in a xy>pea k 
similar to that in a ferromagnetic metal 21J. When k B T raises to a value comparable to 
Wl, further increment in T no longer includes more extended states. Therefore, a xyiPea k 
approaches the universal value of 4 x (ln2)kBe/h. We do, however, observe a smaller a xyiPea k 
value (~18% lower) at N=0(1) near CNP. The broad peak of a xy at N=0(1) is in big 
contrast to the sharp and narrow peak observed in a single-layer graphene at similar field 



strength 



14j . We attribute this to be the lifting of spin degeneracy in the bilayer graphene, 



inferred from the double peak feature in that is attainable at lower field due to higher 
degeneracy at zero-energy LL jg]. A crude estimation of the spin splitting energy AE spin from 
the double-peak in at 15 Tesla gives AE spin ~26 meV (~300 K) which is 15-fold larger 
than the regular Zeeman energy at the same field. In addition, the failure of Mott relation 
near CNP observed in bilayer graphene (Fig. H] (a)) may result from the proposed novel 
phase of counter-propagating edge channels with opposite spin in spin-polarized quantum 
Hall regime [2|. We also remark that the behavior of a xy is well reproduced in S37 and S45 
with sample width to length aspect ratio of 1.3 and 2.4, respectively. Current annealing 
was found to cause insignificant influence except shifting the V g o- The anomaly near CNP, 
arising either from an intrinsic effect of the chiral fermion or from an extrinsic effect such as 



electron-hole puddles 



221 ]. remains an open question. Further investigation in samples with 



higher mobility is required in order to resolve this issue. 

In summary, we performed electric and thermoelectric transport measurements of bilayer 
graphene in quantum Hall regime. Double peak feature in cr xx at N=0(1) LL suggests 
the lifting of the spin degeneracy due to the possible exchange-enhanced Zeeman coupling. 
The disorder width Wl separates two different regimes for the behavior of a xyiPea k- For 
k B T/WL <0.2, a xyiPea k is practically linear in T, giving a xyiPea k/T = 0.19 ± 0.03nA/K 2 
independent of the magnetic field, temperature and LL index N. For k B T/WL >0.5, a xy pea k 
saturates to a value of a xyiPea k — 8.5±1 nA/K close to the predicted value of 4 x (ln2)kBe/h ~ 
9.24 nA/K based on GJ theory of edge current model. We also found anomalous behaviors in 
S yx and a xy near CNP, where semi-classical theory does not provide satisfactory explanation. 
This may imply the existence of a novel phase of counter-propagating edge channels with 
opposite spin in spin-polarized quantum Hall regime, which may have potential application 
in spin-electronics. 
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Note added.- During the preparation of this manuscript, we became aware of a related 
work at lower field (< 7 Tesla) by Nam et al. 23|], which shows consistent result with our 
data. 
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FIG. 1: Thermopower Sxx vs. V g at temperatures ranging from 15 K to 300 K in S37. The 
upper-right inset shows the T dependence of maximum thermopower |S m | and |S xx |(Vg = 80V). 
The red lines are the power law fitting to the data points. The lower-left inset plots 1/S XX vs. V g 
at selected temperatures. 
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FIG. 2: V g dependence of <7 xy and <t xx at temperatures ranging from 8 K to 100 K in S45. An 
optical image of S45 is shown in the upper-right inset. The middle inset plots the sheet resistance 
at CNP Rq vs. T at 15 Tesla (circle) and Tesla (triangle). 
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FIG. 3: V g dependence of (a) S xx and S yx , (b) and <r xy and (c) a xx and a xy at 15 Tesla and 
T=15 K in S37. The dashed lines are the results from S37 after current annealing. The blue thick 
line is a double Gaussian fit to a xx near CNP. 
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FIG. 4: (a) V g dependence of a xy /T (upper panel) and S yx /T (lower panel) in S37 at 15 Tesla and 
three temperatures T=15 K (pentagon), 20 K (diamond) and 25 K (triangle). The orange (upper) 
and green (lower) dashed lines represent the Mott relation fits from d<7 xy /dVg and d(tan#H)/dV g , 
respectively, (b) shows the T dependence of a xy ,peak in S37 (closed symbols) and S45 (open 
symbols) at 15 Tesla and different LLs N=0(1) (circle), 2 (triangle), 3 (star) and 4 (square). The 
lower-right inset plots the a xyiP eak/T vs. ^qH with data points sharing the same symbols as (b). 
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